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We introduce a simple one-parameter network growth algorithm which is able to reproduce a wide 
variety of realistic network structures without having to invoke any global information about node 
degrees such as preferential-attachment probabilities. Scale-free networks arise at the transition 
point between quasi-random and quasi-ordered networks. We provide a detailed formalism which 
accurately describes the entire network range, including this critical point. Our formalism is built 
around a statistical description of the inter-node linkages, as opposed to the single-node degrees, and 
can be applied to any real-world network - in particular, those where node-node degree correlations 
might be important. 



Networks - in particular, large networks with many 
nodes and links - are attracting increasing attention from 
researchers in the fields of physics and biology through 
to sociology informatics and medicine QSSyilSlldH 
Hi IM ^i fill ll^ . Many social and informational net- 
works, including the World Wide Web, exhibit scale- 
free behaviour while others lie closer to fully-random 
or fully-ordered P, Q- Growth mechanisms have been 
proposed for a wide variety of such networks - how- 
ever, many of these growth mechanisms require some 
form of global knowledge of the vertex (i.e. node) de- 
gree. In particular, the famouspreferential-attachment 
model of Barabasi and Albert [j] requires knowledge of 
each node's degree in order to evaluate the corresponding 
attachment probability. Holme and Kim Q and Leary 
et al have subsequently considered a modification of 
this preferential rule by ranking existing nodes according 
to degree, and associating corresponding probabilities of 
connection with the new node accordingly. Their ap- 
proach led to over- and under-skewed distributions for 
the node degree, as compared to the power-law distri- 
bution characteristic of scale- free networks. However in 
both the original and modified preferential-attachment 
algorithms, information is required about the degree - 
and hence preferential-attachment probability - of every 
node in the network, prior to adding a new node or link. 

Although preferential attachment models can yield 
similar distributions to real- world systems, this require- 
ment that the next node or link 'knows' about the degree 
of every node in a large network, makes the mechanism 
microscopically unrealistic for many real-world networks. 
In particular, most biological and social networks are 
too large for such information to be accessible. Moti- 
vated by this shortcoming, the present paper discusses 
a simple, one-parameter algorithm which can reproduce 
under-skewed, over-skewed and scale-free (i.e. power- 
law) networks without global knowledge of the node de- 
grees ~ see Figs. 1 and 2. It therefore provides an alter- 
native, and arguably more realistic, microscopic mecha- 
nism for a wide range of biological and social networks, 
whose growth had previously been explained using the 




FIG. 1: (Color Online) Networks generated by our one- 
parameter, local information growth algorithm with a = 1 
(left), a = 0.2 (center), a = (right). 

Barabasi- Albert preferential-attachment model 0. As 
a by-product of our analysis, we also provide a new for- 
malism to describe general network growth dynamics in- 
cluding node- node linkage correlations. When applied to 
our growth algorithm, this 'link-space' formalism allows 
us to identify the transition point at which over-skewed 
node-degree distributions switch to under-skewed distri- 
butions. It is at this transition point, and only at this 
point, that we find scale- free networks emerge. Simi- 
lar notions of transitive linking have been studied else- 
where - in particular, Saramaki and Kaski and 
Evans |2( have recently presented some fascinating re- 
sults using random-walkers to decide node attachment 
in a network-growth algorithm. However, these analyses 
were all mean-field in nature and hence were not able 
to accurately describe the node-node linkage correlations 
which can be a crucial feature of real-world networks, 
and are a crucial feature of the one-parameter networks 
discussed here. 

We first introduce our general formalism. We then use 
it to discuss various existing growth algorithms and our 
one-parameter growth algorithm, together on the same 
footing. For a growing network in which one node and 
undirected link are added at timestep t, we can write the 
master equation for Xi (the number of nodes of degree 
i) in terms of the probabilities {6j} for attaching to a 
node of degree j : 

\t+i^x, |t+e,_i It -e, It, z> 1 . (1) 

The second term describes the probability of the new 
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Since each link has two ends, the value Lij can increase 
by connection to an (i — l)-degree node which is in turn 
connected to a j-degiee node, or by connection to an 
(j — l)-degree node which is in turn connected to an i- 
degree node. 

We now make a similar steady-state assumption to the 
'node space' example above, and obtain |l3| : 



FIG. 2: (Color Online) Degree distributions for the networks 
in Fig. 1, grown to 100,000 nodes and ensemble-averaged over 
100 networks per a-value. Initial seed comprised two nodes 
and one link. 



node attaching to an existing node of degree i — 1 thereby 
making it a degree i node. The third term on the right 
hand side describes the new node attaching to a node of 
degree i, thereby decreasing Xi. For i — the second 
term is unity as a new node of degree one is added to 
the existing network each timestep. If we assume steady- 
state growth, then Xi = Cit and hence the fraction of 
nodes of degree i is constant — Ci). To facilitate the 
comparison with later expressions, we write the solution 
in the following form: 



1 +^ + ^ 

Cl Cj 



hj > 1; 
J > 1; (4) 



where denotes the fraction of links that connect an 
«-degree node to a j-degree node. It turns out that is 
zero for all the growth algorithms discussed in this paper, 
since the networks that they generate comprise a single 
component. The notation |t has been dropped as before, 
to indicate the steady state. The fraction of i-degree 
nodes is: 
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The notation '|(' has been dropped to indicate steady- 
state. So far we have said nothing about the attachment 
mechanism, and have made the easily geneneralizable re- 
striction that only one node and undirected link is being 
added per timestep. We now follow a similar analysis, 
but retain the node-node linkage correlations that are in- 
herent in many real- world systems 0, ^3 ■ Consider any 
link in a general network - we can describe it by the de- 
grees of the two nodes that it connects. Hence we can 
construct a matrix L such that the element Lij describes 
the number of links from nodes of degree i to nodes of 
degree j. For undirected networks, the summation over 
all elements Lij is equal to twice the total number of The exact solution for k 
links in the network. This matrix represents a surface 
describing the first-order correlations between the node 
degrees - we refer to this as the link space. 

In order to write the master equations for the evo- 
lution of the network in this link space, we must first 
understand the impact of adding a new link. The prob- 
ability of selecting any node of degree i — 1 is given by 
the attachment kernel Suppose an ? — 1 node is se- 

lected - the fraction of these that are connected to nodes 
of degree j, is -jj^jjj^^^^- The term describing the in- 
crease in links from nodes of degree i to nodes of degree 
j, through the attachment of the new node to a node of 
degree i — 1, is given by: 



Hence the degree distribution is retrievable from the link- 
space matrix. To illustrate use of the formalism, consider 
first a random-attachment model in which the existing 
node to which the new node is to be connected, is cho- 
sen randomly. The attachment probability is Qi = Ci. 
Substituting into Eq. [S] we obtain the recurrence rela- 
tion Ci-j_i — Y which yields the familiar 2~' distribution 
of node degree. Substituting into the link-space master 
equation (Eq. Q yields the recurrence relations: 
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1,J > 1, 



(7) 
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where C is the conventional combinatorial 'choose' 
tion. The surface generated is shown in Fig. O 

In the Barabasi-Albert preferential attachment algo- 
rithm 13], the attachment probability is proportional to 
the degree of the node in the existing network: 



X,, 



(3) 
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which in the steady-state hmit of large N (i.e. many 
nodes) can be well approximated by 9^ « Substitut- 
ing this expression into Eq. |21 yields the familiar recur- 
rence relation 



(i - 1) Cj_i 
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whose solution is: 



iii + l)ii + 2) 



(9) 



(10) 



Using the same substitution, the link-space master equa- 
tions yield the recurrence relations: 



ij > 1; 



{i - + (j - l)^»,j-i 

+ i> 1,(11) 

3 + J 



The exact solution for is obtained by algebraic ma- 
nipulation of the link-space matrix using the previously 
derived degree distribution, and is given by [iSj: 

= %W2)r^{G{^ + 1) + 2G(z) - 3G(^ - 1) 
+ 5 El=i(« - l)(a + 6) [G(^ - a) - G(z - a - 1)] }; 



I for a; > 

for X < . 



The first two rows of this matrix, have the form: 



(12) 



2j 



2(6 + .7)(j-l) 
j(j + l)(j + 2)(j+3)' 

2j(j-l)(j + 10) + 48 
3j-(j + l)(j + 2)(j+3)(j+4) 



(13) 



The above equations for lij , even when invoked to low or- 
der in the iteration scheme, can accurately reproduce the 
Barabasi-Albert preferential attachment network. The 
corresponding surface is shown in Fig. |21 

We now turn to our one-parameter network growth 
model, which involves attaching a single node at each 
timestep without prior knowledge of the existing net- 
work structure. The algorithm goes as follows: i) Pick 
a node n within the existing network at random, ii) 
With probability a make a link to that node. Other- 
wise iii) pick any of the neighbours of n at random, and 
link to that node. Hence this algorithm resembles an 
object or 'agent' making a short random-walk. Figure 1 
shows examples of the resulting networks, with Fig. 2 
showing the corresponding degree distributions. Inter- 
estingly, a — Q yields a graph that is dominated by hubs 
and spokes (i.e. extreme over-skewed) 0| while a = 1 
yields the random-attachment graph. Intermediate val- 
ues of a yield networks which are neither too ordered 
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FIG. 3; (Color Online) Representation of the link-space ma- 
trix with elements /ij , in logarithmic coordinates. The val- 
ues shown are the steady-state solutions for the random at- 
tachment (curved surface) and preferential attachment (flat- 
ter surface) algorithms. The value is zero for networks 
consisting of a single component. 



nor too disordered. For a ^ 0.2, the algorithm gener- 
ates networks whose degree-distribution closely resem- 
bles Barabasi-Albert preferential-attachment networks 
(see Figs. 1, 2 and 4). 

We now develop master equations for the evolution of 
this one-parameter network generated with only local in- 
formation. We first establish the attachment probability 
kernel Qi for this algorithm, which in turn requires prop- 
erly resolving the one-step random walk. The link-space 
formalism provides us with an expression for the prob- 
ability P/ associated with performing a random walk of 
length one and arriving at a node of degree i . Note that 
this is different to arriving at a specific node of degree i 
after a one-step walk, since here we consider the possibil- 
ity of arriving at any of the nodes which happen to have 
degree i at time t 



P' = 
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(14) 



This can also be written 0] as P/ = where 
the average is performed over the neighbours of nodes 
with degree i. Note that this quantity does not replicate 
preferential attachment, in contrast to what is commonly 
thought Hill- Defining /3, as 



J2k 



2fe 



yields 



6i = aci + (1 - a)(3iici 



(15) 



(16) 
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FIG. 4: (Color Online) Degree distributions generated using 
the link-space analysis of our one-step algorithm. Inset: re- 
sults for a = 0.25 reproduce the same distribution as for the 
Barabasi-Albert 01 preferential-attachment algorithm. 



Substituting Eq. ^| into Eq. |21 we obtain the following 
for the steady-state node degree: 



,_i(a+(l-a)/3,_i(z~l)) 
1 a + (1 - a)|3^i 



(17) 



Substituting into Eq. 0] yields the following recurrence 
relations: 



+ (1 - a)/3i_i(i - 1) 



+ (1 - a)/3,-_i(i - 1) 



1 + 2a + (1 - a)(i3i + jfij) 



1,3-1 



+ (1 - 0)^3-1 (j - 1) 



+ (1 - a)flj„iO - 1) 



1 +2o + (1 - a)(,3i 

The non-linear terms resulting from [3 imply that an ex- 
plicit closed expression for lij is difficult. We leave this 
as a future challenge - but we stress that our formalism 
can be implemented in its non-stationary form numeri- 
cally (i.e. iteratively) with very good efficiency 0, as 
demonstrated by the degree distributions in Fig. 4. 

Our numerical results suggests that for a ~ 0.2, the de- 
gree distribution approaches scale-free. We now use the 
above formalism to deduce the critical value Oc at which 
the node-degree distribution goes from over-skewed to 
under-skewed, and hence the value of a at which a scale- 
free distribution arises. At Oc, we know that the attach- 
ment kernels for both our and the Barabasi-Albert model 
should be equal. Hence from Eq. |Sland Eq. ^1 we have 



(18) 



which for large i yields ac — 1 — We could proceed 
to use the exact solution of the link-space equations for 



the preferential-attachment algorithm, in order to infer 
/3i in the high i limit. However, since can be expanded 
in terms of li^^ as shown in Eq. El and j decays very 
rapidly as i, j become large, we can obtain a good ap- 
proximation by only using the first two terms of Eq. 1131 
This gives /3 ~ 0.66. Substituting into Eq.^jthen yields 
the critical value at which scale-free networks exist as 
flc = 0.25, in excellent agreement with the results shown 
in the inset plot of Fig. 4. 



In conclusion, we have presented a simple, one- 
parameter algorithm for generating networks whose prop- 
erties span from the exponential degree distribution of 
random attachment, through to an ordered 'hub and 
spoke' situation (Fig. 1). A scale-free network turns out 
to be a special case, lying at a transition point in between 
the two limits. The growth algorithm utilizes only one 
simple parameter a, and requires no global information 
concerning node degree. As a by-product of the analysis 
of this model, we have managed to develop a new type 
of link-space formalism which can account for the node- 
node linkage correlations in real- world networks. Indeed, 
as will be shown elsewhere, the link-space formalism can 
be used to describe an even wider variety of networks 
than those emerging from the present growth model. 
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For mathematical details and derivations, see 
www.physics.ox.ac.uk/users/smithdmd/hnkspace.pdf 
Suppose we seed the algorithm with an existing hub and 
spoke graph, and use a = 0. The probability that a new 
node will be linked to the hub tends to unity Eis the net- 
work grows. More generally, the distribution of vertex 
degrees will approach that of a hub and spoke network 
as the network grows. This type of network can be con- 
sidered the most ordered in that the mean free path is 
minimized (i.e. it tends to two) 



